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We observe a signal for the doubly charmed baryon H^ c in the charged decay mode 
H^ c — > K~ir + in data from SELEX, the charm hadro-production experiment at Fermilab. We 
observe an excess of 15.9 events over an expected background of 6.1 ± 0.5 events, a statistical sig- 
nificance of 6.3a. The observed mass of this state is 3519 ± 1 MeV/c 2 . The Gaussian mass width of 
this state is 3 MeV/c 2 , consistent with resolution; its lifetime is less than 33 fs at 90% confidence. 

PACS numbers: 14.20.Lq, 13.30. Eg 
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The addition of the charmed quark to the (uds) triplet 
extends the flavor symmetry of the baryon octet and 
decuplet from SU(3) to SU(4). Even though the large 
m c breaks the symmetry, SU(4) still provides a good 
classification scheme for baryons composed of uds and 
c pj. There is strong experimental evidence for all the 
predicted baryon states which contain zero or one va- 
lence charmed quark jlj. In this letter we present the 
first experimental evidence for one of the six predicted 
baryon states which contain two valence charmed quarks 
- the doubly charmed baryons. There have been many 
predictions of the masses and other properties of these 
states @-|f|. The properties of doubly charmed baryons 
provide a new window into the structure of baryonic mat- 
ter. 

The SELEX experiment uses the Fermilab 600 GeV/c 
charged hyperon beam to produce charm particles in a 
set of thin foil targets of Cu or diamond. The three- 
stage magnetic spectrometer is shown elsewhere ||[§. 
The most important features are: (a) the high-precision, 
highly redundant, silicon vertex detector that provides 
an average proper time resolution of 20 fs for single- 
charm particle decays, (b) a 10 m long Ring-Imaging 
Cherenkov (RICH) detector that separates ir from K up 
to 165 GeV/c ||, and (c) a high-resolution tracking sys- 
tem that has momentum resolution of ap/P < 1% for a 
200GeV/c reconstructed A+. 

The experiment selected charm candidate events us- 
ing an online secondary vertex algorithm. A scintillator 
trigger demanded an inelastic collision with at least four 
charged tracks in the interaction scintillators and at least 
two hits in a positive particle hodoscope after the second 
analyzing magnet. Event selection in the online filter re- 
quired full track reconstruction for measured fast tracks 
(p>15 GeV/c). These tracks were extrapolated back into 
the vertex silicon planes and linked to silicon hits. The 
beam track was measured in upstream silicon detectors. 
A full three-dimensional vertex fit was then performed. 
An event was written to tape if all the fast tracks in 
the event were inconsistent with having come from a sin- 
gle primary vertex. This filter passed 1/8 of all interac- 
tion triggers and had about 50% efficiency for otherwise 
accepted charm decays. The experiment recorded data 
from 15.2 x 10 9 inelastic interactions and wrote 1 x 10 9 
events to tape using both positive and negative beams. 
67% of events were induced by E~, 13% by 7r~, and 18% 
by protons. 

This analysis began with a sample of A+ single- 
charm baryons decaying to pK~ir + . Candidates were 
selected with a topological identification of 3-prong 
positively-charged secondary vertices, requiring a mo- 
mentum measurement for each track. RICH identifica- 
tion of the proton and kaon was required. Charged tracks 
with reconstructed momenta which traversed the RICH 
[p > 22 GeV/c) were identified as protons or kaons if those 
hypotheses were more likely than the pion hypothesis. 



The other positive track was identified by the RICH as 
a pion when possible, otherwise it was assumed to be a 
pion. The primary vertex was refit using all other found 
tracks. Details of the selection procedure can be found 
in [ |io||ri| ]. This A+ sample (1630 events) contains the 
same A+ — ► pK~ir + events we used to measure the A+ 
lifetime pi. 



A Cabibbo-allowed decay of a doubly charmed baryon 
must have a net positive charge and contain a charmed 
quark, a strange quark and a baryon. We chose to search 
for decay modes like — » A+ K~tt + with an intermedi- 
ate K~ir + secondary vertex between the primary vertex 
and the A+ vertex. In this analysis we have incorpo- 
rated the charm selection techniques developed for single- 
charm baryon states, as above. 
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FIG. 1. A candidate event with production target, la track 
error corridors and vertex error ellipses. This is a plan view of 
three dimensional tracks and vertices. Three additional found 
tracks which form the primary vertex with the beam track are 
not shown. 

Events were analyzed for evidence of a secondary ver- 
tex composed of an opposite-signed pair between the 
primary and the A+ decay point. We used all tracks 
not assigned to the A+ candidate in the search. A new 
primary vertex was formed from the beam track and 
tracks assigned to neither the A+ nor the K~tt + ver- 
tices. The new secondary vertex had to have an accept- 
able fit x 2 arL d a separation of at least la from the new 
primary. The A+AT~7r + transverse momentum with re- 
spect to the incident beam direction is required to be 
in the range 0.2 < p t [GeV/c] < 2.0. These cuts were 
developed and fixed in previous searches for short-lived 
single-charm baryon states. We have applied them here 
without change. 

Most tracks from the K~tt + vertex have insufficient 
momentum to reach the RICH. Masses were assigned 
according to topology. For the signal channel negative 
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tracks are assigned the kaon mass and positive tracks the 
pion mass. As a background check we also kept wrong- 
sign combinations in which the mass assignments are re- 
versed. A candidate event from the A+A' _ 7r + sample is 
shown in Fig. |l|. Further details of the A+ reanalysis may 
be found in Ref. [Hi. 
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FIG. 2. (a) The A+K~n + mass distribution in 5 MeV/c 2 
bins. The shaded region 3.400-3.640 GeV/c 2 contains the 
signal peak and is shown in more detail in (c). (b) The 
wrong-sign combination A^K + 7v~ mass distribution in 5 
MeV/c 2 bins, (c) The signal (shaded) region (22 events) and 
sideband mass regions with 162 total events in 2.5 MeV/c 2 
bins. The fit is a Gaussian plus linear background. 

In Fig. ||(a) we plot the invariant mass of the A+A~7r + 
system, fixing the A+ mass at 2284.9 MeV/c 2 Q. The 
data, plotted in 5 MeV/c 2 bins, show a large, narrow 
excess at 3520 MeV/c 2 . This excess is stable for different 
bin widths and bin centers. Fig. |](b) shows the wrong- 
sign invariant mass distribution of the A+A + 7r~ system 
with the same binning as in (a). There is no significant 
excess. 

In Fig. ||(c) the shaded region from (a) is re-plotted 
in 2.5 MeV/c 2 bins and fit with a maximum likelihood 



technique to a Gaussian plus linear background. The 
fit has x 2 /dof = 0.45, indicating that the background is 
linear in this region. 

To determine the combinatoric background under the 
signal peak we exploit the linearity of the background jus- 
tified by the fit. We define symmetric regions of the mass 
plot in Fig. |(c): (i) the signal region (3520 ± 5MeV/c 2 ) 
with 22 events; and (ii) 115 MeV/c 2 sideband regions 
above and below the signal region, containing 162 — 22 = 
140 events. We estimate the number of expected back- 
ground events in the signal region from the sidebands as 
140 * 5/(115) = 6.1 ± 0.5 events. This determination has 
a (Gaussian) statistical uncertainty, solely from counting 
statistics. The single-bin significance of this signal is the 
excess in the signal region divided by the total uncer- 
tainty in the background estimate: 15.9/^(6.1 + 0.5 2 ) = 
6.3cr 0. The Poisson probability of observing at least 
this excess, including the Gaussian uncertainty in the 
background, is 1.0 x 10~ 6 . 

Our reconstruction mass window is 3.2-4.3 GeV/c 2 
with 110 bins of width 10 MeV/c 2 in this interval. The 
overall probability of observing an excess at least as large 
as the one we see anywhere in the search interval is 
1.1 x 10~ 4 . 

This state has a fit mass of 3519 ± 1 MeV/c 2 . Our 
expected mass resolution, from a simulation of the decay 
E+ -> A+A~7r+ is ~ 5 MeV/c 2 . We observe a Gaus- 
sian width of 3 ± 1 MeV / c 2 , consistent with our simula- 
tion. The confidence level for a fit with a Gaussian width 
fixed at our expected resolution is 20%. The width we 
observe is consistent with statistical fluctuations in this 
small sample. 

The wrong-sign mass combination is plotted in 
Fig. ||(b). Those events show comparable fluctuations 
to the sidebands of the signal channel and give no evi- 
dence for a significant narrow structure. We have inves- 
tigated all possible permutations of mass assignments for 
the non-A+ tracks. The peak at 3520 MeV/c 2 disappears 
for any other mass choice, and no other significant struc- 
tures are observed. Reconstructions with events from the 
A+ mass sidebands produce relatively few entries and no 
significant structures in the doubly charmed baryon re- 
gion. 

The dependence of the signal significance, as defined 
above, on several of the selection cuts is shown in Fig. |^. 
The significance depends strongly only on the A~7r + ver- 
tex separation. The dependence is driven by a large in- 
crease in background at small separations and the ab- 
sence of both signal and background events at large sep- 
arations. The dependence on the width of the signal 
region is stable, only decreasing when made wider than 
the mass resolution. All other cuts have been checked; 
no significant dependence on any cut has been seen. 

A weakly-decaying state has two c quark decay 
amplitudes plus a W-exchange amplitude for c + d — > 
s + u. This suggests that its lifetime will be of the order 
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FIG. 3. The signal significance as a function of: (a) the 
vertex significance of the the K~tt + vertex, (b) the half width 
of the signal region, (c) the mass window around the A J mass, 
(d) the half width of the sideband region. The square points 
are the values used in this analysis. 

of the S° or shorter, rather than like the long-lived S+. 
For each event we compute the reduced proper time 
t = (m/p) x (L — L m i n )/c. Here L is the flight path from 
the primary vertex to the Sj c decay point, and L m j„ 
is the minimum flight path for a given event. In this 
analysis L min = a, the error on the flight path. One can 
form the likelihood, £(r) = Yli^i P(*i; r ) to observe a set 
{t{\ of reduced proper time values, given the probability 
P(t,r) for a given lifetime r. P(t,r) is the normalized 
probability of observing a given reduced proper time t for 
a true lifetime r determined from the SELEX simulation. 
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TABLE I. The change in likelihood, —log(C) + log(£ max ) 
for different Ej. lifetimes. Increases of 0.5 and 2.0 in this 
statistic correspond to la and 2a respectively. 

Results for the change in likelihood from its maximum 
value are given in Table [j] for both the signal region and a 
sideband region, where combinatorial backgrounds dom- 
inate. In both cases the lifetime is much shorter than the 
single-charm S° lifetime of about 100 fs Q. Given our 
single-event resolution of 20 fs |l(|, we cannot exclude a 
prompt decay for the signal events from the likelihood. 
However, the fact that the significance of the signal in- 



creases as the flight path L is increased [see Fig ||(a)] 
suggests that the lifetime is not zero. From these data 
the upper limit for the S+ lifetime at 90% confidence is 
33 fs. 
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TABLE II. Events produced by incident particle species. 
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fraction of AJ signal events 


0.68 


0.32 
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sideband region events 
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47 





TABLE III. Events produced by different target materials. 

Production characteristics of the 22 signal plus back- 
ground events are indistinguishable from the single- 
charm A+ sample. |ll| The mean p t is 1 GeV/c and 
mean xf ~ 0.33. It is interesting to compare produc- 
tion of the S+ state by different beam hadrons. Results 
for the signal region and sidebands are listed in Table 0. 
One sees that the doubly charmed baryon candidates are 
produced solely by the baryon beams. The top row gives 
the fraction of all interactions produced from each of the 
beams. The proton/E _ ratio for production and for 
all interactions is the same. One can also ask if there a 
dependence on the target nucleus. 



Table III shows that 
the diamond/copper ratio of the signal events is similar 
to the sideband events, which in turn behave like single- 
charm production. Production of candidates differs 
from A+ production only in the absence of meson beam 
contributions. 

The yield of this state is larger than most production 
models predict jl2|]. The acceptance for the 15.9 events 
we observe in this final state, given that we observe a A J, 
is 11%. Using a factor 1.5 from isospin to account for the 
S+ — > A+K 7T° mode and Bjorken's estimate [[l3| of 1.6 
to include other decay modes with A+ in the final state, 
we find that ^20% of the A+ in this sample are produced 
by H+ decay. 

Our production region has not been probed by other 
experiments. The large xp, small p t region is not 
amenable to perturbative QCD analysis. Two CERN 
experiments with tt~ beams have reported anomalously 
large production of events with two charmed parti- 
cles [ ^4|Jl~5| ] . However, we do not know of any model calcu- 
lation that would predict this large hadro-production rate 
of doubly charmed events. The FOCUS charm photo- 
production experiment at Fermilab has searched for dou- 
bly charmed baryons in their charm samples. They see 
very few candidate events and no signal peaks ES]. The 
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BELLE B-factory experiment recently reported a very 
large production ratio for (J/ipcc)/(J/tp) in continuum 
e + e~ annihilation at 10.8 GeV p7[ . The connection to 
the hadro-production data is not clear. 

In summary, we have observed a narrow state at 3520 
MeV/c 2 decaying into A+K~ir + , consistent with the 
weak decay of the doubly charmed baryon 2+ . We report 
this state as the first observation of a doubly charmed 
baryon. 
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